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Abstract 

Neutralinos that are mostly wino or higgsino are shown to be compatible with the 
recent DAMA annual modulation signal. The nucleon scattering rates for these dark 
matter candidates are typically an order of magnitude above the oft-considered bino. 
Although thermal evolution of higgsino and wino number densities in the early universe 
implies that they are not viable dark matter candidates, non-thermal sources, such as 
from gravitino or moduli decay in anomaly mediated supersymmetry breaking, suggest 
that they can be the dominant source of cold dark matter. Their stealthiness at high 
energy colliders gives even more impetus to analyze nucleon scattering detection meth- 
ods. We also present calculations for their predicted scattering rate with germanium 
detectors, which have yet to see evidence of WIMP scattering. 
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1 Introduction 



There are two intriguing reasons to study the lightest neutralino from supersymmetric ex- 
tensions of the Standard Model. For one, we would like an answer for what composes the 
non-baryonic matter necessary to explain the observed rotations of galaxies. By angular 
momentum considerations, we note that the observed visible matter in our galaxy implies a 
lower theoretical angular velocity than observed. We also observe galaxies in clusters insuf- 
ficiently bound by the visible matter. Secondly, we demonstrate in this article that should 
light neutralinos exist as a cold dark matter (CDM) candidate, the detection prospects are 
promising for a wide range of parameter space, including for higgsinos and winos. 

In the various supersymmetric models, the sneutrino, the gravitino, and a neutralino 
are all capable of being the lightest supersymmetric particle (LSP) Q. Here, we study 
the lightest neutralino \ as the LSP. Depending on the supersymmetry breaking model 
considered, the lightest neutralino may be bino-like, wino-like, higgsino-like, or admixtures 
of these three categories. 

Historically, the bino-like LSP has been studied extensively as a dark matter candidate |2|. 
The wino-like and higgsino-like LSP have been largely ignored as dark matter candidates 
because they can be demonstrated in the majority of popular models, that their thermal 
relic densities are simply too low for dark matter contention ||. 

We include the possibility of a relatively new class of models — the anomaly-mediated 
supersymmetry breaking (AMSB) scenarios @, |l| ||. The relevant characteristic of these 
models is the possibility of a wino-like or higgsino-like LSP (from a non-thermal source) 
with a density capable of being the dark matter |7], [8]. The wino-like LSP is not unique to 
the AMSB scenarios. Theoretical grounds are also provided by moduli-dominated supersym- 
metry breaking within O-II superstring motivated models 13] and supersymmetry breaking 



in which the responsible F-term is not a SU(5) singlet [ I0| . What is unique about AMSB is 
its prediction that the gravitino mass is a few orders of magnitude above the superpartner 
masses. This characteristic is what ultimately enables the wino or higgsino to be a viable 
dark matter candidate. 

If the lightest neutralino truly is the dark matter or a significant fraction, they, being 
weakly interacting in character, will penetrate the earth much like a neutrino. This allows 
for possible detection by scattering elastically with nucleons. Of the different neutralino 
LSP detection methods |fLl|| , we report on the detection prospects of the lightest neutralino 
scattering off nucleons within different target samples. 
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2 Super symmetry Breaking Models 



First, a brief review of how the neutralino obtains its character is supplied. The Standard 
Model includes an electrically neutral U(l) gauge field B and a neutral SU(2) gauge field 
W 3 . The minimal supersymmetric extension includes two neutral Higgs scalar fields, H® 
and H®. These four fields have supersymmetric fermionic partners if) = (B, W 3 , H%, H®) 
with various interactions amongst one another, allowing for a non-diagonal mass matrix to 
be formed out of the Lagrangian quadratic terms. A transformation may be performed to 
diagonalize the neutralino mass matrix to the physical basis, x — {xu X2, X3, X'di with xi 
having the smallest mass eigenvalue. The Haber and Kane convention QT2Q for fi is used in 
the mass matrix. The notation x will De used to mean the lightest neutralino If, for 

example, the largest component of x is the W 3 component, we say the LSP is wino-like. 
Thus, the LSP's eigenvector components determine its character, as made explicit by 



X = N n B + N 12 W 3 + N 13 H d + N l4 H u . (1) 

We consider gravity mediated and anomaly mediated supersymmetry breaking. The val- 
ues of M\ and M 2 are determined within these models, with /1 constrained only by additional 
assumptions on the superpartner spectrum. In a large class of minimal supergravity super- 
symmetry breaking scenarios, scalar masses mo and gaugino masses M, are chosen to be the 
same at the unification scale. In these models, the relation between the U(l) and SU(2) 
gaugino masses is 



— y = — y ( m i/2 universal at GUT scale). (2) 

This amounts to M\ = | tan 2 Qy/M<i being about half M 2 at the electroweak scale. Within 
these models, there is still a chance that \x is smaller than M\. The lightest neutralino x will 
then be either bino- or higgsino-like for these models. 

A wino-like x is possible within the AMSB class of models. AMSB models supply two 
results relevant to our discussion. The first of which claims that the gaugino masses are 
proportional to the gauge coupling beta function, 

Mn = QtLm 3/2 (AMSB) (3) 

9n 

where g n is the gauge coupling constant, m 3 / 2 is the gravitino mass, and (3 gn is the beta 
function of the gauge coupling constant. At the electroweak scale, the values of g n and (3 9n 
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result in M2 being approximately a third of M\. Furthermore, \i < M2 is possible, depending 
on the sparticle spectrum. In this case, AMSB would produce a higgsino LSP. 

The wino's cross-section for pair annihilation is higher than the bino's for two reasons: 
g 2 is larger than g 1; and bino pair annihilation has less options for the final state than winos. 
For example, winos can annihilate to W pairs without any intermediate scalar superpartners 
whereas pure binos cannot. This results in a present day number density riy, too low for winos 
as a cold dark matter candidate (but often just right for binos, which originates its appeal as 
a cold dark matter candidate). This rules out a wino-like LSP and is the reason the wino-like 
LSP has been largely ignored as a CDM candidate. Let us call this a "thermal source" for 
wino production. This disucssion also applies to an AMSB-produced higgsino-like LSP. 

AMSB supplies a second, "non-thermal source" for \ production. The gravitino or mod- 
uli fields in this model are what constitute the second LSP source via direct or indirect \ 
decay products. A modulus field may be defined as the scalar field that parameterizes a flat 
direction in the theory. The moduli fields may acquire a mass from supersymmetry breaking 
as the moduli's own potentials are "lifted" at the ends. The important feature of AMSB mod- 
els is the granting of a large mass (10 to 100 TeV) to the gravitino and moduli fields. These 
large masses can circumvent the historically annoying gravitino problem and "cosmological 
moduli problem" in which the gravitino and moduli acquire small masses under common 
supergravity models. A small enough mass may have a lifetime longer than successful Big 
Bang nucleosynthesis would allow. In general, it is undesirable for any cosmological scenario 
to include a field dumping Plank scale energy into the universe after nucleosynthesis begins. 
The large masses from AMSB models ensure the gravitino and moduli decay sufficiently 
early 0,| . 

3 The Scattering Rate 



The detection prospects of a neutralino from the galactic halo are quantified as the elastic 
scattering (or event) rate R. The assumed virialized neutralinos are bound to the galactic 
halo by a gravitational potential, and the scattering rate (measured in events/day/kg) is [|T3|] 

R = ap * V ft. (4) 

where = m^n^ is the mass density of the neutralinos, is the average speed of the 
neutralinos as they float around the galactic halo, Mjy is the mass of the target nucleus, and 
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(a) h- and H -exchange (b) Z-exchange 

Figure 1: The leading order parton diagrams for iV — \ scattering 



c) ^-exchange (d) scattering with gluons 
via heavy quark loops 



F% is the nuclear form factor (we take F^ — 1). In terms of the assumed average values, the 
neutralino's speed and density, the scattering rate is 



R 



rrtyM N 



^-1.8 x 10 n GeV 4 



Pic 



events 



v 0.3GeV/cmV V320km/sec/ \day-kg, 



(5) 



The remainder of this section is organized as follows: our model's lagrangian is stated, the 
possible processes are discussed, and the many coupling constants are made explicit. The 
Lagrangian is divided into two parts: spin-dependent terms and spin-independent terms. 
This separation is motivated by realizing spin-independent interactions are enhanced by the 
presence of all the nucleons in the nucleus. 

To describe the spin-independent interactions, we use an effective Lagrangian interaction 
term of the form 

£C fN[M][NN), (6) 

Ne{p,n} 

where the subscript N is the nucleon field, proton or neutron, is the scalar 4-point effective 
coupling constant that includes Higgs and squark exchange, as well as the neutralino-gluon 
scattering of Fig. 1(d). The contribution of scattering with gluons occurs through a loop 



of quarks or squarks and is fully calculated in Ref. JL4J]. It is shown to be at least an 
order of magnitude smaller than Higgs exchange in Ref. |13] for most of parameter space 
involving reasonably heavy squarks. For this reason, only one neutralino-gluon interaction 
is large enough for our task — Higgs exchange to a triangle of heavy quarks (c, b, t) coupled 
to gluons. These interactions are in addition to another effective Lagrangian term that 
describes spin-dependent Z, Higgs, and squark exchange: 



(7) 



where c q and d q are effective coupling constants to be made explicit shortly. 
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These Lagrangian terms yield a scattering cross-section of [TJj 



a = — 



m x M N 
my + Mj 



N 



(n p f p + n n f n f + 4X 2 J(J + 1) ( £ d q Aq 



The factor A 2 J(J+ 1) represents the fraction of the nucleons' spin in comparison to the total 
nucleus spin (squared) J(J + 1), which includes orbital angular momentum. 



and 



n n are the masses of the lightest neutralino, a proton, and a neutron, respectively. Aq is a 
quark's second moment of the quark density of its spin polarization |T5| . 

The spin-independent effective coupling constant fx is composed of two terms: 



f N = fu + f D - (9) 

The spin-independent squark exchange contributions are organized into terms proportional 
to the sums and differences of MSSM Yukawa coupling constants a qi and b qi from the q~iqx 
vertex, supplied in the Appendix. The terms proportional to a qi — b qi are represented by 
fu- Those proportional to a qi + % are small and omitted /# is the spin-independent 



contribution due to Higgs exchange to all six flavors of quarks. The form of /# is 

f H 2 f H 

fu= £ —fT q m N + — £ —fTgrn N (10) 

q€{u,d,s} m i Z ' q£{c,b,t} m i 

where fog is the fraction of the nucleon mass the light quarks (u, d, s) effectively represent, 
defined by fTq m N = (N\m q qq\N) for nucleon N. fo g = 1 — 2 fr q is the remaining fraction 
carried by the nucleon sea. f q is the effective coupling constant for a /^[xx][9<z] Lagrangian 
term that describes Higgs exchange, f q follows from a straightforward MSSM calculation 
of h and H t-channel exchange between x an d a quark: 



fH _ m j C hxx C hqq _|_ CH X X C Hqq \ 

q Q \ rn 2 h m\ J 



(11) 



where the MSSM Yukawa coupling constants for Higgs-neutralino and Higgs-quark are sup- 
plied in the Appendix. 

The squark exchange interactions proportional to a~. — 6? have an effective coupling of 

fo= £ ^f Tq m N . (12) 

q&{u,d,s} m i 
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This refers to simple s-, t-, and -u-channel exchanges with the light quarks. All squark 
exchange interactions proportional to a~ —b]- , denoted fs, are at least an order of magnitude 
smaller than the total interaction /, even for reasonably light squarks. 

Squark exchange appears as both spin-dependent and spin-independent. The 4-point 
Lagrangian interaction term implied by s-, t-, and -u-channel diagrams of neutralino-quark 
scattering by squark exchange is [qx] [xq] . The Fierz rearrangement supplies the scalar form 

ccfumm (13) 

as well as the vector /axial- vector form 

>Cc[| 7 M 7 5 x][^(4 + 47 5 )?], (14) 
where /|, c|, and S q are effective coupling constants: 



1 2 a? -6? 
Jq " 4^m|-(rn* + m 9 ) 2 ^ > 



1 dnb 



d = 



y ■^g 1 / 16 x 

2^ 1 mi-(m* + m q y 

and b qi refer to the MSSM couplings for the qqx vertices and are supplied in the Appendix. 

As for spin-dependent interactions, Z-exchange and squark-exchange involve a vector 
coupling constant c q and an axial vector coupling constant d q : 

c q = 4 + c z q: d q = dl + d z q . (18) 

c z and d z are effective coupling constants for a 4-point xxqq vertex that implicitly proceeds 
through Z-exchange: 

9 2 2 OUT q 3-2e q sm 2 e w ) z _ gjO'^ ( . 

c q ~ /i™2 ' a q ~ /i_2 v iy / 



where T q s is the quark's weak isospin third component, e q is the electric charge of the quark, 
and 0" R is the Zxx vertex coupling from the MSSM (see the Appendix). 
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Figure 2: Contributions to R (events/day /kg) vs. m x (GeV). These are the contributions to 
the x _73 Ge scattering rate R by Higgs (solid), squark (dotted), and Z exchange (dashed). 
These contributions are calculated from (a) a mSUGRA bino-scenario (Mi = M 2 = §/■*) 
, (b) an AMSB wino-scenario (Mi = M 2 = and (c) a higgsino-scenario (Mi = |//, 
M2 = 3/i). The remaining parameters are tan/3 = 4, squark masses rrtq = 2 TeV, and soft 
trilinear parameters Aq = 0, pseudo-scalar Higgs mass uia = 500 GeV, and the nuclear and 
astrophysical parameters of Table 1. 

4 Scattering off Germanium 



We now numerically compare the contributions to the scattering rate by the Higgs, squark 
and Z exchange in Figures 2. The target sample is germanium to match the CDMS 



and HDMS |T7| collaborations. The nuclear and astrophysical parameters used are listed in 
Table 1. The nuclear parameters are from the 1988 EMC data M, |H| |20|] and Ref. [2~I|. For 
reasonable parameter constraints supplied by AMSB and mSUGRA, we created scenarios 
with the lightest neutralino as wino-like, bino-like, and higgsino-like. From the numerical 
results, the Higgs-exchange contribution is seen to be much larger than squark-exchange 
and Z-exchange. The contribution from Z-exchange becomes significant for low values of 
\i (which is roughly proportional to 171%) due to the cross-section contribution of the form 
m r j /i 4 where m r is the reduced mass of the nucleus and neutralino. 

To analytically demonstrate the Higgs dominance, first consider the effective coupling 
constant of the 4-point interaction that proceeds through Higgs exchange. It is pro- 
portional to the MSSM h\X an d H\X vertices. In the approximations of wino-like lightest 
neutralino {N12 ^> Nu) and small Higgs mixing angle a, these coupling constants are 
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c hxx w ~92N 12 N 14 (20) 

CH XX » -^2^12^13- (21) 

Further imposing the limits mz/\fJ> ± M\\ 1, mz/\fJi ± M 2 | 1, and tan/? > 2, these 
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Table 1: Nuclear and astrophysical parameters used for calculations [18, 1£, 20, pjfl . 



for protons: 




fau 


0.023 


fad 


0.034 


fas 


0.14 


for neutrons: 




fau 


0.019 


fad 


0.041 


fas 


0.14 



Au 


0.77 


Ad 


-0.49 


As 


-0.15 


for 73 Ge, A 2 J(J + 1) 


0.065 




1 




320 km/s 




0.3 GeV/cm 3 



coupling constants become 



g 2 M 2 + fj, sin 2/3 
c hxx - y /i 2 — Mf mw (22) 



# 2 // cos 2/3 

w y^Mf mw - (23) 



For sketching purposes, these couplings go as ~m vl ///i. The contribution of Higgs-exchange 
to the cross section will then go as ~m^m 2 v //i 2 m^. 

The coupling, , for xq scattering through Z-exchange (Eqn. [H]) goes as ~m|//x 4 and 
contribution to the cross-section goes as ~m 2 m^//i 8 . The wino- and higgsino-like neutralino 
scenarios of figures 2(b) and 2(c) show over an order of magnitude Z-exchange increase over 
the bino-like scenario since their cross-sections are proportional to g%, versus g\. 

The squark-exchange couplings /|, c|, and d^ are proportional to a~. + a^b^, and 
a \ ~ ^li - a ?i an< ^ ^ De °^ ^ e same order for all types of neutralinos if tan/3 is not 
extremely large or small. An extreme value of tan (3 may cause squark-exchange to become 
more pronounced (see the Appendix for the Z q0 coupling). For moderate values of tan/3, 
it is because the squark cross-section contribution is proportional to ~ rn%/rri~ that squark 
exchange is suppresed when a heavy squark is considered. 

Fig. 3 shows the numerically calculated scattering rates for neutralinos and 73 Ge at 
tan/3 = 4 for the three cases of "moderate" values of the squark mass (m, = 2 TeV) and 
pseudo-scalar Higgs masses (tua = 500 GeV), a light squark (rrig = 500 GeV) case, and a 
light pseudo-scalar Higgs (rriA = 150 GeV) case. The practical aspect of these plots lies in 
recognizing the various possibilities of tuning the cross-section and scattering rate by varying 
the supersymmetry breaking model parameters /x, m^, and rriA along with tan/3. The rates 
for 76 Ge are similar to those for 73 Ge since the spin-dependent contributions are small for 



8 



0.1 1 




0.1 



0.01 



0.001 



0.1 



0.01 



0.001 



100 200 300 400 



(a) £>-like x 



100 200 300 400 



(b) W-like x 



100 200 300 400 



(c) H-like x 



Figure 3: The x- 73 Ge scattering rates R (events/day/kg) as a function of the LSP mass m x 
(GeV) for bino, wino, and higgsino scenarios. Scenario (a) uses M\ = Af 2 = §/U, (b) 
uses Mi = M 2 = and (c) uses Mi = M 2 = 3fi. The solid lines have "moderate" 
parameters rriq = 2 TeV, rriA = 500 GeV. The dashed lines represent a light squark with 
nig = 500 GeV, = 500 GeV. The dotted lines represent a light pseudo-scalar Higgs with 
rriq = 2 TeV, = 150 GeV. For all, tan/3 = 4 and soft trilinear parameters Ag = are 
used. 



both. 

The bino-like neutralino scattering rate suffers from the lower g\ coupling constant, ex- 
plicit in Eqs. ^ and 27. The wino- and higgsino-like neutralinos have comparable scattering 
rates due to the common vertices that originate from the Lagrangian terms — ^(^iJjifjH^ 3 . 
The bino-like neutralino has analogous terms but are proportional to the smaller U(l) gauge 



coupling g\. Analogous to Eqs. |22] and £3|, using the same limits, the higgsino-like LSP from 
mSUGRA has the Yukawa couplings 



g 2 1 + sin 2(3 , . 

°hxx ~ — a — m w tan%, (24) 

4 fi — Nl\ 



9i 2 sin 2 (5 
4 (j, - Mi 



c Hxx ~ — ^7 — m w tanv w . (25) 



While these calculations show promising event rates for a wide range of parameter space, 
the regions of parameter space that result in undetectable scattering rates should be noted. 
Larger values of tan/3 will lower the hit coupling (Eq. ^) and suppress the light Higgs 
exchange. Large tan/3 also slightly increases the light Higgs mass and causes further 
suppression. In the case of light squarks, a large tan f3 has the opposite effect — enhancing 
the scattering rate by causing significant left-right squark mixing, driving the sbottom mass 
down. The corresponding squark propagators are then enhanced. Overriding the effects of 
varying tan /3, larger mass parameters for the Higgs, gauginos, and squarks can quickly reduce 
the scattering rate to undetectable levels. When tan/3 3> 1 the cross-section begins to rise 
again. One can show that the cross-section starts to rise like tan 2 /3 when tan/3 > m H/ m ^ 
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Figure 4: Illustrations of the DAMA constraints (roughly equivalent to the boxes in each 
figure panel). Each plot represents the scalar (spin- independent) contribution to the cross- 
section (Tscaiar (pb) as a function of the wino-like neutralino mass (GeV) for neutralino- 
proton scattering in an AMSB scenario (Mi = M 2 = |/i). (a) represents "moderate" 
parameters (rriq = 2 TeV, ttia = 500 GeV), (b) represents a light squark scenario (m, = 500 
GeV, rriA = 500 GeV), and (c) represents a light pseudo-scalar Higgs scenario (rriq = 2 TeV, 
rriA = 200 GeV). A q = is used for all soft trilinear parameters. The values of tan/3 are 
2 (solid), 10 (dashed) and 30 (dotted). The upper (lower) set of plots have values // > 

Gu<o). 

due to the increased relative importance of the heavier Higgs boson. 



5 DAMA Constraints on Super symmetry Parameters 



The DAMA collaboration reports indications of a WIMP with a mass of 52^ u GeV and 
a cross-section of £c" sca i ar = 7-2-o.i X 10~ 6 pb when scattered against a single proton, using 
standard astrophysical assumptions (as in Table 1). £ is the fraction of the cold dark matter 
that this WIMP represents. The "scalar" subscript on the cross-section refers to only the 
spin-independent contributions. 

We apply our analysis to interpret this WIMP signal as a scattering of the lightest 
neutralino. While the CDMS collaboration claims to exclude the neutralino parameter space 



implied by DAMA, the DAMA limits are used here at face value [IB|. The AMSB wino 
scenario is chosen as the framework for this interpretation without too much loss of generality. 
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One can roughly equate this with the higgsino scenario, or take off an order of magnitude 
from the scattering rate to imagine the mSUGRA bino scenario (as in Fig. 3). 

We have made assumptions about the astrophysical and squark parameters (Ag and nig). 
Altering the astrophysical models, i.e. the galactic halo structure, may significantly change 



the allowed SUSY parameter space. Refs. [23, 24, 25 1 demonstrate a means for different 
galactic halo models to allow for WIMP masses up to 150 GeV at 3a. Non-universal squark 
masses and soft trilinear parameters allow for a broad range of cross-sections for any fixed 



set of values m^, tan/3, and sgn(/i), as seen in Refs. |26], |27], |28|, P9fl . As for assumptions 



on the nuclear parameters, Ref. |3(| demonstrates the effect of uncertainties of the quark 
masses and mass fractions fx q - 

Using the "moderate" SUSY parameters of Fig. 4(a) {rriq = 2 TeV, rriA = 500 GeV, 
Ag = 0) the cross section is significantly below the favored region of the DAMA/NaI-2 data. 
The lighter squark of Fig. 4(b) will increase the scalar cross-section, but not to DAMA 
heights. A light pseudo-scalar Higgs is necessary for compatibility with the DAMA limits. 
By tuning tan /?, rriq, and rriA the DAMA interpretation may be made compatible with 
anomaly and gravity mediated SUSY breaking models. 

One concern with the above parameter space with light pseudo-scalar is that the lightest 
CP-even Higgs mass will fall below the LEP-II limits. Indeed, our lightest Higgs mass in 
our examples range from about 100 GeV to 115 GeV which is just at the LEP-II SM Higgs 
limits. However, one should keep in mind that the LEP-II SUSY Higgs limits are significantly 
below the SM Higgs limits due to decreased hZZ coupling when the pseudo-scalar is light. 
Furthermore, our result for matching the DAMA allowed region is mostly independent of 
the light Higgs mass or the squark masses. Therefore, the rest of the superpartner masses 
can arrange themselves to produce large enough loop corrections to satisfy the Higgs mass 
constraints without significantly affecting our final result. 



6 Conclusions 

We complement previous efforts to demonstrate the compatibility of supersymmetry and the 
recent DAMA annual modulation signal [p7|-[35| by showing that the wino-like LSP from 



AMSB models are also compatible with the DAMA signal. Furthermore, the AMSB models 
lead to detectable event rates for a large volume of SUSY parameter space through the tuning 
of tan/3, rriq, and rriA- The AMSB-inspired models produce nearly identical event rates if 
a higgsino is the LSP. This is due to the dominance of Higgs exchange and the Hi «-> W 
symmetry in the —^g2HiHiW 3 operators. 
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mSUGRA and AMSB models differ by about an order of magnitude in the event rate 
when mSUGRA produces a bino-like LSP, due mainly to the differences of the SU(2) and 
U(l) gauge couplings. The wino-like neutralino's event rate is always higher than the bino- 
like neutralino's for all of parameter space. 

Although the WIMP mass and cross-section implied by DAMA's recent observations are 
confirmed to be compatible with mSUGRA and AMSB models, the necessary parameters are 
near the edge of what is ruled out by other experiments, such as collider physics. For wino 
and higgsino LSPs, the masses can be rather light (certainly lower than 80 GeV) without 
running into collider constraints because it is so difficult to find winos and higgsinos at 
colliders 131 
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Appendix: MSSM Coupling Constants 

For the Higgs-fermion-fermion vertices, using the notation of Drees and Nojiri (Ref 

Chxx = ^(#2^2 - giN n ) (sin aiVi 3 + cos aN u ) (26) 

chxx = 2^ 92Nl2 ~ 9i N n)(- cos aN n + sinaA^i 4 ) (27) 
92 sin a 

Chdd - „ {2o) 

2mw cos p 

92 cos a 

c H dd = —z n 29 

2mw cos p 
92 cos a 

Chun = —Z — T (30) 

2mw sm p 
92 sin a 

Chuu — ~- : — n \ 61 ) 

2m w sm p 

where a is the Higgs mixing angle. 

For the Z\X vertex: 



r 2 



\Nu\ 2 - \N 13 \ 2 



(32) 



For the qqx vertices: 



an = -[cos9g(X q0 + Z qQ ) +sin9g(Y q0 + Zgo)) (33) 
a<j 2 = g [ _ sin 0q( x <io + Z qo) + c °s 9q(Y q0 + Z q0 )] (34) 
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bqi = ^[cos^X^q - Z q0 ) + sm9q(-Y q0 + Z q0 )} (35) 
bq 2 = -[- sin 6g(X q0 - Z q0 ) + cos ^(-Kgo + z q o)} (36) 
where #g is the squark mixing angle of left and right squarks into physical squarks and 

X q0 = -V2g 2 [T g3 N l2 -Un9 w (T q3 -e q )N n } (37) 
Y q0 = V2g 2 %tme w e q N xl (38) 
i 2m "f 14 for up-type quarks 

g2in d h 
\/2 cos /3m w 
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